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Abstract 

Multiplicity fluctuations in limited segments of momentum space are calculated for a classical pion gas 
within the statistical model. Results for the grand canonical, canonical, and micro-canonical ensemble 

o , 

\ are obtained, compared and discussed. We demonstrate that even in the large volume limit correlations 
"q ■ between macroscopic subsystems due to energy and momentum conservation persist. Based on the micro- 
canonical formulation we make qualitative predictions for the rapidity and transverse momentum dependence 
CN \ of multiplicity fluctuations. The resulting effects are of similar magnitude as the predicted enhancement due 
' ' to a phase transition from a quark-gluon plasma to a hadron gas phase, or due to the critical point of strongly 
' interacting matter, and qualitatively agree with recently published preliminary multiplicity fluctuation data 
^ of the NA49 SPS experiment. 
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I. INTRODUCTION 



The statistical model has been, for a long time, successfully applied to fit experimental data on 
mean hadron multiplicities in heavy ion collision experiments over a wide range of beam energies 
and system sizes. For recent reviews see . So naturally the question arises whether the 

statistical model is able to describe event-by-event fluctuations of these observables as well. And 
indeed, a first comparison suggests that this might be possible for the sample of most central events. 
Global conservation laws, imposed on a statistical system, lead, even in the large volume limit, to 
suppressed fluctuations. The multiplicity distributions of charged hadrons recently reported [5| by 
the NA49 SPS experiment are systematically narrower than a Poissonian reference distribution. This 
could be interpreted gI as effects due to energy and charge conservation in a relativistic hadronic 
gas. 

Multiplicity fluctuations are usually quantifled by the ratio of the variance of a multiplicity dis- 
tribution to its mean value, the so-called scaled variance. In statistical models there is a qualitative 
difference in the properties of mean value and scaled variance. In the case of the mean multiplicity 
results obtained within the grand canonical ensemble (GCE), canonical ensemble (CE), and micro- 
canonical ensemble (MCE) approach each other in the large volume refers here to as the 
thermodynamic equivalence of these ensembles. It was recently found [7] that corresponding results 
for the scaled variance are different in different ensembles, and thus this observable is sensitive to 
conservation laws obeyed by a statistical system. 

The growing interest in the experimental and theoretical study of fluctuations in strong interac- 
tions (see e.g., reviews Sj) is motivated by expectations of anomalies in the vicinity of the onset 
of deconflnement and in the case when the expanding system goes through the transition line 
between a quark-g.uou plasma and a hadron gas phase R In particular, a cr.tical point of strongW 
interacting matter may be accompanied by a characteristic power-law pattern in fluctuations 
A non-monotonic dependence of event-by-event fluctuations on system size and/or center of mass 
energy in heavy ion collisions would therefore give valuable insight into the phase diagram of strongly 

interacting matter. Provided the signal survives the subsequent evolution and hadronization of the 
fl 

system (see also uM)- Therefore, in order to asses the discriminating power of proposed measures, 

n n n n n 

for a recent review see pj[, one should flrstly study properties of equilibrated sources [6|, lU, 115|, 116[ 
and quantify 'baseline' (or thermal/statistical) fluctuations. Apart from being an important tool in 
an effort to study a possible critical behavior, the study of fluctuations within the statistical model 
constitutes also a further test of its validity. 

In this paper we make detailed predictions for the momentum space dependence of multiplicity 
fluctuations. We show that energy and momentum conservation lead to a non-trivial dependence of 
the scaled variance on the location and magnitude of the observed fraction of momentum space. These 



predictions can be tested against existing and future data from the heavy ion colhsion experiments 
at the CERN SPS and BNL RHIC facihties. 

The paper is organized as follows: In section [TTl we briefly introduce our model. In section IIIII we 
consider multiplicity distributions in a limited region of momentum space in GCE and CE. For the 
MCE we follow, in section IIVI the procedure of Ref. 17| and show how to calculate the width of the 
corresponding distributions in the large volume limit. We revisit the so-called 'acceptance scaling' 
previously suggested as an approximate implementation of experimental acceptance in section IVl 
Technical details of the calculations are presented in the Appendix. Concluding remarks and a 
summary in sections |Vl] and IVIII close the paper. 



II. THE MODEL 

The ideal Boltzmann 7r+ vr^ vr*^ gas serves as the standard example throughout this paper, while the 
main subject of investigation is the multiplicity distribution P{Nq) of particles with momenta inside 
a certain segment Q of momentum space. Calculations are done for the three standard ensembles 
GCE, CE, and MCE. For the sake of argument we will assume that we only want to measure P{N^), 
i.e. the probability distribution of negatively charged pions in a limited segment Q of momentum 
space. Hence vr" with momenta inside Q are observed, while n~ inside the complementary segment 
fl are not observed. 7r+ and vr^ are never detected. In GCE and CE the presence of 7r° as a degree 
of freedom is of no relevance, while in MCE it constitutes a heat bath for the remaining system. For 
consistency we use the same system throughout this discussion. 

In order to keep the model simple, we assume a static homogenous fireball. Our considerations 
therefore exclude collective motion, i.e. flow, and resulting momentum spectra are purely thermal. 
We also omit resonance decay contributions in this work. The spectra presented in Fig. [1] are 
normalized to the total vr" yield in GCE and CE. Thus they are the same in both ensembles. In 
MCE one expects in the large volume limit only small deviations from Boltzmann spectra. None 
of the forthcoming arguments are affected by this. In the following we will use the transverse 
momentum and rapidity spectra presented in Fig. [T]to construct bins fli = ApTi = [pt,,Pt,+i] (left), 
or Qi = Ai/i = [yi,yi+i] (right), as indicates by the drop-lines. 

In section UTT] we calculate the multiplicity distributions P{N^) for arbitrary segments Q for the 
ideal Boltzmann GCE and CE. To characterize the distribution one can calculate its (raw) moments 
{N^) from: 

oo 

(N^) = E NnP{Nn) • (1) 

Nn=0 
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FIG. 1: Differential particle spectra for a classical pion gas at T = 160MeV. Left: Transverse momentum 
spectrum Eq. ijBlh . Right: Rapidity spectrum Eq. ijCip . Both curves are normalized to the total yield Eq.([6]). 
The bins are constructed in a way that each bin contains 1/9 of the total yield. 

A convenient measure for the width of a distribution is the scaled variance: 



(Nn) 



(2) 



In order to remove simple scaling effects, the bin sizes or segments are chosen such that each bin 
or segment contains the same fraction q = (Nq) / {N^t^) of the total yield (compare Eq.([2])). Here 
(Nq) denotes the average particle number in the momentum space segment Q, and (A^47r) denotes 
the average total {4tt integrated) multiplicity. The effect of finite acceptance can approximately be 
taken into account by 



UJn 



1 + q {Uj4n - 1) 



(3) 



where 004^.,^ assumes the ideal situation when all particles are detected, while Uq assumes that particles 
are detected with probability q regardless of their momentum. Hence Eq.Q holds when particles 
are assumed to be uncorrelated in momentum space. In the limit g — > one observes a random 
distribution with Ug 1, i.e. a Poissonian, while when q 1 one sees the real distribution with 
width ujq uo^T^. In this work we take explicitely correlations due to globally conserved charge (CE), 
and energy- momentum (MCE) into account and compare the results to Eq.([3]). 



III. GRAND CANONICAL AND CANONICAL ENSEMBLE 



A. Grand Canonical Ensemble 



In the GCE, both, heat and charge bath are assumed to be infinite. And thus neither charge, en- 
ergy nor momentum are conserved exactly. Temperature T and charge chemical potential /i regulate 
average energy and charge density in a system of volume V . Usually it is said that charge, energy 
and momentum are conserved in the average sense and fluctuations about an equilibrium value are 



allowed. Apart form Bose and Fermi effects 18| particles are therefore uncorrelated in momentum 
space. However this example serves as an illustration for the following CE and MCE calculations. 
We start by decomposing the Boltzmann single particle partition function z~ {(pNn) of ^~ i^^^o 
parts, 



^~ (</>iVo) = (M,) = J d'pe — e'^-o + J d'pe " , (4) 

where the single particle energy e = ^/p'^ + m^, and m, and g are mass and degeneracy factor of 
7r~ respectively. Only for momentum states inside the momentum space region Q we introduce 
additionally a Wick-rotated fugacity exp {itp^^)- For the positive and neutral pion (which we do not 
want to detect in our example) we write: 



^ , d p e T and z = ^ / d p e t . (5) 

i2nY J i2nf J ^ ^ 

The value of the single particle partition function, for instance of the neutral pion, is given by: 

z' = {N') = ^^m'TKJ^). (6) 



27r VT 

For the sake of simplicity we assume equal masses for all pions. To obtain the GCE multiplicity 
distribution for in a momentum space segment Q we use the Fourier integral over the generalized 
GCE partition function Z {(Pnq) = exp [z^ {(pNn) ~^ ^ci ~^ ~^ ^^~\ ' non^alized by the GCE partition 
function: 

P,,ANn) - Z^^lx |^e-^"^-oZ(0^J = M^I^^ exp , (7) 

— TT 

where the system partition function is given by Zg^e = Z {(t>NQ = 0), and = {4>Nn = 0). Inde- 
pendent of the shape or size of Q we find a Poissonian for the multiplicity distribution Eq.Q. Thus, 
using Eq.(l2]), one finds for the scaled variance ujf^^ = 1, since {Nq) = z^, and (A*"^) = {Nn)^ + {N^). 

For Bose and Fermi statistics one does not expect a Poisson distribution and (in particular when 
the chemical potential is large) deviations from a Poissonian can be large. Thus one expects also 
deviations from Eq.Q when considering only finite acceptance. 



B. Canonical Ensemble 



In the CE the heat bath is still assumed to be infinite, while we remove the charge bath and drop 
the chemical potential. Thus, we introduce a further Wick-rotated fugacity fi/T ^ i(f)Q into the single 
particle partition functions to account for global (however not in the momentum space segment VL) 
conservation of electric charge Q. Particles in Vt are therefore correlated, due to the condition of fixed 



net-charge, with a finite charge bath composed of 7r+ and unobserved vr . We again spht the single 
particle partition function for tt~ into an observed, (^atj^, </)q), and an unobserved part, 



gv 



while we do not want to measure tt"*" and 7r°, and thus: 

gV 



{2n) 



dp e Te 



and 
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gV 



d^p e 
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d^p e" 



(9) 



The normalization of the CE multiplicity distribution is given by the CE system partition function 
Zcei i.e. the number of all micro states with fixed charge Q, Z^^ = Jg (2z) exp(z°), where Iq is 
the modified Bessel function. The multiplicity distribution of Nq in a momentum space segment Q, 
while charge Q is globally conserved, can be obtained from Fourier integration of the generalized GCE 
partition function Z {(pM^^, 0q) = exp [z^ {(pj^^, 0q) +2;^ (0q) + z+ (0q) + 2°] , over both angles 4>q 
and 



( -\Nn 00 / xa 
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where in CE Zr 



{<pNn =4>Q = 0), % = 



{Q + Nn + a 
0), and z = 2;+ 



respective first two moments one finds from Eq. 

- Iq+1 i2z) 



and 



{Nl) = {z^) 
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Thus, we obtain the well known canonical suppression of yields |19l . 
The result, however, is completely independent of the position of the segment Q. And therefore the 
scaled variance, Eq.([2]), takes the form: 



l + z; 
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Iq (2^) 



(13) 



where uq is the width of Pce{NQ), i.e. the multiplicity distribution of vr" with momenta inside Q, 
while u^T^ is the width of the corresponding distribution when Q is extended to the full momentum 
space. It can immediately be seen that this formula is consistent with acceptance scaling, Eq.(l3]), 
Ufi = 1 + q {uj^T, — 1), if g = z^/z. Generally we find ujI'^. < uj^ < uj^'^^ = 1. In the limit of z^/z 
we approach the Poisson limit of a 'random' distribution with u = 1, i.e. the observed part of the 
system is embedded into a much larger charge bath and the GCE is a valid description. 



IV. MICRO-CANONICAL ENSEMBLE 



For the MCE an analytical solution seems to be out of reach presently, so we use instead the 
asymptotic solution, applicable to large systems, derived in [17|]. In order to avoid unnecessary 
repetition of calculations, we will only give a general outline here, and refer the reader for a detailed 
discussion to Ref-jisj. It should be mentioned that this method be would be also applicable to 
systems of finite spatial extention, provided the average particle number in a given momentum space 
bin exceeds roughly {Nq) > 5. In this work we confine ourselves to large systems and try to asses 
the general trends. 

The basic idea is to define the MCE multiplicity distribution in terms of a joint GCE distri- 
bution of multiplicity, charge, energy, momentum, etc. The MCE multiplicity distribution is then 
given by the (normalized) conditional probability in the GCE to find a number Nq of particles in 
a segment Q of momentum space, while electric charge Q, energy E, and three momentum P are 
fixed. Therefore we will keep temperature and chemical potentials as parameters to describe our 
system. Effective temperature and effective chemical potential, i.e. Lagrange multipliers, can be 
determined by demanding that the GCE partition function is maximized for a certain equilibrium 
state (Q, E, P). This requirement is entirely consistent |l7] with the usual textbook definitions of T 
and /i in MCE and CE through differentiation of entropy and Helmholtz free energy with respect to 



conserved quantities. In principle we would have to treat all conservation laws on equal footing |24l |. 
and thus introduce Lagrange multipliers for momentum conservation as well. However here we are 
only interested in a static source, thus P = 0, and the relevant parameters are equal to zero. 

In the large volume limit energy, charge, and particle density in the MCE will correspond to GCE 
values. This is required by the thermodynamic equivalence of ensembles for mean quantities. MCE 
and CE partition functions are generally obtained from their GCE counterpart by multiplication with 
delta-functions, which pick out a set of micro states consistent with a particular conservation law. 
Here it will be of considerable advantage to use Fourier representations of delta-functions, similar to 
the treatment in Section IIIII This method could be considered to be a Fourier spectral analysis of 
the generalized GCE partition function [itI ]. 

The normalized conditional probability distribution of multiplicity Nq can be defined by the ratio 
of the values of two partition functions: 

_ number of all states with Nq, Q, E, and P = , , 

Pmce[Nn) = — — — — . (14) 

number of all states with Q, E, and P = 
The real MCE partition function and our modified version are connected as Z{V, Nq,Q, E, P) = 
Z^"''^'^'^(V,T, ^)e~^^^^e~^'^^^ . In either case the normalization in Eq.( fT4|) is given by the partition 
functions with fixed values of Q, E, P, but arbitrary particle number A^^^, hence Z{V, Q, E, P) = 
E^„=o Z{V, Nn, Q, E, P), or Z«'^'^(\/, T, /i) = ^^^^^ ^A.o,Q,i?,P(^^ j^^ However when taking the 



ratio ([13]) auxiliary parameters chemical potential and temperature drop out: 

^ Z{V,Nn,Q,E,P) ^ Z^-^<^'^'^{V,T,fi) 

Z{V,Q,E,P) ZQ'^'P{V,T,f,) 

The main difference between the two versions of partition functions is that for Z{V,Nq,Q,E,P) 
one is confronted with a heavily oscillating (or even irregular) integrant, while for Z^^^'^'^'^iV^T^ ^) 
the integrant becomes (T,/i correctly chosen) very smooth. Thus, introduction of T and /i allows to 



derive (and use) the asymptotic solution of Ref . [Jj 

We have a total number of 6 conserved 'charges', and hence we need to solve the 6-dimensional 
Fourier integral for the numerator in Eg. (1151) ^: 
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X exp 



(16) 



y^i'k (0iVo0Q, 0p,, 0Pj,, 0pj 

k 

The summation in f|T6|) should be taken over the single particle partitions Vipk = Zk of all considered 
particle species k. The Wick-rotated fugacities 0q, etc. are related to the individual conservation 
laws. The distinction between the Kronecker delta-function (limits of integration [— tt, tt]) for discrete 
quantities and the Dirac delta-function (limits of integration [—00,00]) for continuous quantities is 
important here, however for deriving an asymptotic solution it will not be. To simplify flTB]) we change 
to shorthand notation for (pj = {(pN^-^cpg, (j)E, (pp) and the conserved 'charge' vector = {N^, Q, E, P). 
We again split the single particle partition functions in two parts. The first part counts the number 
of momentum states observable to our detector, while the second part counts momentum states 
invisible to our detector: 

(0,) = ^JSp e^«i"<^^ + d'p e"^ e<^'^^ . (17) 

(27r) J (27r) J 

For the 'charge' vector of all measured particle species k we write ^1^ = {l^qk-,Sk-,f>k) for momenta 
inside and (fi.Q = (0, g^, £^fc,Pfc) for momenta outside of Q. For all unobserved particle species we 
write Q ~ lin~ (0' Iky ^k,Vk)- Here qu is the electrical charge of particle species fc, and Sk and pk 
are its energy and momentum vector. In Ref. 17|, where only multiplicity distributions in the full 
momentum space were considered, the general 'charge' vector took the form g;^^^ = {nk,qk,£k,Pk), 
where is the multiplicity of this particle. For stable particles = 1 in case they are observed, 
and nfc = if they are not measured, while for unstable particles could also denote the number 
of measurable decay products. 



We drop in the following the argument {V, T, fi) to simplify the notation. 



For large system volume the main contribution to the integral (fT6|) comes from a small region 
around the origin [isj]. Thus we proceed by Taylor expansion of the integrant of f|T6l) around (f)j = 
0. In this context (0j) = J^k'^kifpj) would be called the cumulant generating function (CGF). 
Cumulants (expansion terms) are defined by differentiation of the CGF at the origin: 



= -2 



;i8) 



<t>j=0 



Generally are cumulants tensors of rank n and dimension defined by the number of conserved quan- 
tities. Here ki is a 6 component vector, while ^2 is a 6 x 6 matrix, etc. 

The parts of the integrant related to discrete quantities, i.e. and Q, are now not 27r periodic 
anymore (while in Eq. f|T6|l they are), but superpositions of oscillating and decaying parts. Thus we 
extent the limits of integration to ±oo, what introduces a negligible error. Eq. (|T6l) therefore simplifies 
to: 



n 

i=i. 



exp 



n=0 



n 



(19) 



Summation over repeated indices is implied. Existence and finiteness of the first three cumulants 
provided, any such integral can be shown to converge to a multivariate normal distribution in the 
large volume limit: 

exp 



'gee 



(20) 



'(27ry)''/^det |a| ' 

where Zg^e = exp [Vkq] is the GCE partition function, is the cumulant of O*'* order, = 
[Q^ — V {^^^)k y^^^"^ is a measure for the distance of a particular macro state to the peak 
Vk^I of the joint distribution, and a is the square root of the second rank tensor ^2, see [171] for 
details. 

The normalization in Eq. flT^ can essentially be found in two ways. The first way would be to 
integrate the distribution (|20|) over all possible values of multiplicity N^, while all other variables 
are set to their peak values, e.g. Q = Vk^ , E = Vnf, P = 6. The second and more practical 
way is to use an approximation similar to Eq.( !20l) to describe the macro state = {Q, E, P). The 
normalization in Eq.( JT5l) . Z^'^'^ , is then given by the 5-dimensional integral, similar to Eq.( JT6l) . 
without the integration over ipN^- The 1-dimensional slice along N^, i.e. the conditional distribution 
of particle number A^^f^, while charge, energy and momentum are fixed to Q,E,P = 0, can then be 
shown [17] to converge to a Gaussian in the large volume limit: 

(iVQ-(iVn))'' 



Pmee{Nf, 
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(27r (Nn)) 



1/2 



exp 



o , ,rnee 



(21) 



The scaled variance uJq'^'^ is given by the ratio of the two determinants of the two relevant second 
rank cumulants, k,2 and k2, of the two partition functions Z^^^'^'^'^ and Z^'^'^ , hence^: 

det 1^2 1 



UJr 



K 



det \k2\ 



(22) 



The asymptotic {V —>■ oo) scaled variance can therefore be written in the form of Eq.(28) in 17l |. 
Considering only the asymptotic solution we need to investigate only the first two cumulants (n = 
1,2) in detail. We will first discuss the structure of ni and ^2, and then deduce a few properties of 

The first order cumulant ki of z'^^'^'^'^ gives GCE expectation values for particle density , 
charge density , energy density nf , and expectation values of momentum k^''", etc. Since we are 
only interested in a static source we find due to the antisymmetric momentum integral (see Appendix 



0. The general form of the first cumulant Ki is then: 



Hi = Kf, U, U, Uj . (23) 

The second cumulant K2 of Z^'^'^'^'^ contains information about correlations due to different 
conserved quantities. A detailed discussion of correlation terms only involving Abelian charges 
and/or energy, e.g. k^'*^, k^'^, and k^'^, can be found in jl?!. Again, due to the antisymmetric 
nature of the momentum integral, all cumulant entries involving an odd order in one of the momenta, 
or K^'^"" are equal to zero. The general second order cumulant ^2 thus reads: 
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(24) 



Please note that by construction, Eq. fllSI) . the matrix (l24l) is symmetric, hence k^"' 



K, 



etc. 



The second matrix k2, now related to the partition function Z^'^'^, is obtained from K2, Eq.(l2l 
by crossing out the first row and first column. In the following we are going to make use of the fact 
that one can express the determinant of a matrix A by: 

n 

det|A| = 5^(-l)^+'A,,fcM,-fc , (25) 

where Aj^k is the matrix element j, k of a general non-singular n x n matrix A, and Mj^k is its 
complementary minor. A simple consequence of Eq. fl25l) is: 

2" 



det \ k2\ 
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Px,Px ,Py,Pv ,^Pz,Pz 



K 
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E,E Q,Q 
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f/€^-'^-)^det|/€2| 



(26) 



^ Please note, that in order to simplify formulas, the notation is slightly different from [T 



where n^^'^" = ^2^'^^ = k^'"^" > spherical symmetry in momentum space, and /t2 is just a 2 x 2 

matrix involving only terms containing E and Q. In case correlations between particle number and 
conserved momenta are vanishing, i.e. k^'*"'^"' = 0, or k^^'^"" = 0, then, similarly to Eq. (!26|) . the 
determinant of k,2 factorizes into a product of correlation terms {k^'"^'')'^ the determinant of a 
3x3 sub-matrix involving only terms containing E, Q, and A^. Hence in taking the ratio Eg. (1221) one 
notes, that in this case momentum conservation will not affect multiplicity fluctuations in the large 
volume limit jl^]. In this work, however we do not necessarily find k^"'^"" = 0, as we only integrate 
over a limited segment Q of momentum space, and taking momentum conservation into account may 
affect the result. 

Finally it should be stressed that this procedure can be easily generalized to account for Bose or 



ermi statistics. Also phenomenological phase space suppression (enhancement) factors 7g |25|] or 7, 



261] could be straightforwardly included. However, without proper implementation of the effect of 
additional correlations due to resonance decay and collective motion, i.e. flow, it seems of little value 
to do too strict calculations for experimentally measurable distributions. We thus return to the pion 
gas example from section IIIII and restrict the discussion to simple momentum space cuts in rapidity, 
transverse momentum, and azimuthal angle, see also the Appendix for details. 

V. RESULTS 

A. Multiplicity fluctuations in the full momentum space 

Let us firstly recall basic properties of multiplicity fluctuations of negative particles in the full 
momentum space {4tt fluctuations) in the three standard ensembles, of the Boltzmann pion gas 
considered here. 

Multiplicity fluctuations in the CE are suppressed due to exact charge conservation. For a neutral 
[Q = 0) system one finds in the large volume limit u^'^ = 0.5 7|. Further suppression of fluctuations 
arise from additionally enforcing exact energy conservation in the MCE. Here one finds u^^'^ ~ 0.25 
for a Boltzmann pion gas at T ^ IGOMeV. In the GCE, since no conservation laws are enforced, we 
always find a Poisson distribution with width ool^^ = 1. 

Since charge conservation in CE links the distributions of negatively charged particles to the one 
of their positive counterparts, i.e. P{N_) = P{N^ — Q), the relative width of P{N_) increases 
(decreases) as we move the electric charge density to positive (negative) values [27!]. This can be 
easily be seen from Eq. flMl) by crossing out all rows and columns containing energy and momentum 
and calculating the asymptotic scaled variance of negatively charged particles, u^, from Eq. (!22i) . 



The same effect is present in the MCE, however the calculation is slightly longer. 

Results for Arc multiplicity fluctuations of negatively charged particles in a Boltzmann pion gas at 
T = 160MeV and different charge densities are summarized in Table [H Additionally estimates, based 
on our previously employed 'uncorrelated particle' approach, Eq(l3]), for multiplicity fluctuations with 
limited acceptance are given. Despite the fact that u^-j,- is very different in GCE, CE, or MCE and 
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TABLE I: Multiplicity fluctuation of vr^ in a classical pion gas in the large volume limit in the three standard 
ensembles at T = 160MeV for different charge densities. The index '47r' denotes fluctuations in the full 
momentum space, while the index 'g = 1/9' assumes acceptance scaling, Eq.(l3]). The ratio n^/ritot equals 
to 0.33 for /X = 0, 0.48 for /x = -m/2, and 0.20 for ^ = +m/2. 

also rather sensitive to the charge density, the estimates for limited acceptance (g = 1/9) based on 
Eq.([3]) vary only by a few %. In order to decisively distinguish predictions for different ensembles a 
large value of q would be needed. 

B. Multiplicity fluctuations in limited segments of momentum space 

In Section IIIII we have seen that in the Boltzmann CE multiplicity fluctuations observed in a 
limited segment of phase space are insensitive to the position of this segment. The dependence on 
the size of the segment can thus be taken into account by use of acceptance scaling Eq.(l3]). To 
balance charge a particle can be produced or annihilated anywhere in momentum space. And due to 
a infinitely large heat and momentum bath in the CE no momentum state is essentially preferred. 

In the MCE this dependence is qualitatively different. When using the MCE formulation particles 
are correlated due to the constraints of exactly conserved energy and momentum, even in the large 
volume limit. Fluctuations in a macroscopic subsystem are strongly affected by correlations with the 
remainder of the system. 

In Fig. [2] we show the scaled variance of multiplicity fluctuations for negatively charges particles 
in finite bins in transverse momentum (left), and rapidity (right). The bins are constructed such 
that each bin contains on average the same fraction q of the total average yield. The width of each 
bin is indicated by the bars. Calculations are done for two values of acceptance (g = 1/5, and 
q = 1/9). The dashed and dotted lines correspond to acceptance scaling Eq.Q, while the markers 



are calculated from Eq. (!22|) . One finds that multiplicity fluctuations in bins with high transverse 
momentum and high values of rapidity are, due to energy and momentum conservation, essentially 
suppressed with respect to bins where individual particles carry less energy and momentum. 
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FIG. 2: Transverse momentum (left) and rapidity dependence (right) of the scaled variance of vr" at T = 
IGOMeV , for a classical pion gas at zero charge density. Momentum bins are constructed in a way that each 
bin contains the same fraction q of the average 7r~ yield. The (error)bars indicate the width of the pT or y 
bins, while the marker indicate the center of gravity of the corresponding bin. The lines indicate acceptance 
scaling Eq.([3]). Calculations are done for different values of acceptance, q = 1/5 (square marker, dotted 
line), q = 1/9 (triangle down, dashed). 

A intuitive explanation would probably look like this: Let us consider an event with an unusually 
large (small) number of particles at the most forward rapidity bin. In this bin we would find therefore 
a macroscopic state with unusually large (small) observed longitudinal momentum P°'"' and energy 
j^obs^ The remainder of the system therefore has to have rather large (small) momentum —P°^-^ 
and rather small (large) energy E — E"^"^. Since both probability distributions, for the observed 
and the unobserved subsystems, do not factorize into independent probabihty distributions, but are 
correlated, this macro state would be rather unlikely. Fluctuations about the mean (A^^,) at forward 
(backward) rapidities should therefore be suppressed. On the other hand can modest multiplicity 
fluctuations in a high pt bin induce stronger fluctuations in the lower pt bins, and fluctuations about 
(Apj,) in a low px bin are enhanced. Even when detecting only a fraction of about 10% of the total 
system these correlations can have a sizeable effect. 



C. Conservation laws 



It seems worthwhile to consider individual conservation laws and their impact on multiplicity 
fluctuations in more detail. On of the main advantages of the analytical procedure presented here, is 
certainly that one can easily 'switch on' or 'switch off' a particular conservation law. For illustrative 
purposes we show the result of cu^^^ for MCE without longitudinal momentum conservation in FigJS] 



In comparing Figj2l right, to Figj3] it becomes obvious that energy conservation alone cannot 
account for the strong suppression of multiphcity fluctuations at forward rapidities, but has to be 
explained by combined energy and longitudinal momentum conservation. 



without P, conservation 





»q = 1/5 




= 1/9 




FIG. 3: Same as Figl2] right, however without Pz conservation. 



The relevant cumulants elements, which give information about the strenght of correlations be- 
tween particle number and a particular conserved quantity, 

(I'Px etc. Whenever 

a element is vanishig, then the corresponding conservation law has no impact on multiplicity fluc- 
tuations. For details of the calculations please see the Appendix. Since for fluctuations of charged 
particles k^"'*^ and k^"''^ are generally non-zero, we will focus only on the effects of momentum 
conservation. 

For multiplicity fluctuations in bins in transverse momentum momentum conservation does not 
affect the result, see Appendix [HI and the suppression effect is a result of energy conservation alone. 



When considering cuts in rapidity one finds in general k^^'^^ 7^ 0, but K2 



0, and only 



longitudinal momentum conservation needs to be taken into account, see Appendix O In considering 
the third idealized case, where our detector observes only a segment in azimuthal angle 0, but all 
rapidities y and transverse momenta pt, both global P^, and Py conservation lead to non-trivial 
modifications of Eq.(l3]), see Appendix [Dl 

To understand the difference between the strong suppression of fluctuations at high transverse 
momentum and the rather modest suppression at high rapidity when momentum conservation is not 
enforced, one should compare the elements k^*"^'^ in Eq. flB3|) . and k^^''^ in Eq. flC4l) . which measure 
in Boltzmann approximation the average energy density carried by particles in a bin Ap^ or Ay, 
to the total average energy density {E_) = /t^*'"'^ carried by vr^. (All other elements in Eqs. flB3l) 
and (1C4I) do not depend on the location of the segment.) In case of kinematical cuts in ApT the 
fraction k,2^^' /(-^-) rises from about 5% in the lowest to roughly 20% in the highest p^-bin. In 
contrast to that for the central y-bin this ratio is about 10%, while the most forward or backward 
bins it is roughly 12%. However in both cases the bins contain on average g = 1/9 ~ 11% of the 
total average vr" yield. The effect of energy conservation is thus weaker for cuts in rapidity than for 



cut in transverse momentum, see also Appendices [B] and O 



D. Charged systems 

In Figs, mthe transverse momentum (left) and rapidity (right) dependence of the scaled variance 
is presented for two different values of charge density. Similar to the CE, in MCE the effective size 
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FIG. 4: Transverse momentum (left) and rapidity dependence (right) of the scaled variance of vr" at 
T = 160MeV, for a classical pion gas at zero charge density. Momentum bins are constructed in a way that 
each bin contains the same fraction q of the average vr" yield. The (error)bars indicate the width of the 
Pt or y bins, while the marker indicate the center of gravity of the corresponding bin. The lines indicate 
acceptance scaling Eq.(l3]). Calculations are done for different charge densities /x = — y (square marker, 
dotted), and = + y (triangle down, dashed). 



of the heat and charge bath matters. We find that in general MCE effects for negatively charged 
particles are stronger (weaker) when the electric charge density is negative (positive). In the limit 
of a strongly positively charged system, the vr" subsystem could be considered as embedded in a 
large heat, charge, and momentum bath (provided by 7r+ and 7r° particles) and MCE effects would 
cease. The GCE would here be the appropriate limit. In the opposite limit of a strongly negatively 
charged system, charge conservation essentially becomes equivalent to particle number conservation. 
This scenario might be more familiar from textbooks, where the CE is usually understood as the 
ensemble with fixed particle number. However here also the same arguments as above apply, except 
the effect would be much stronger, and u^^^ = 0. 

In general one would expect that suppression effects in bins of high transverse momentum or high 
values of rapidity are stronger the more abundant the analyzed particle species is. In the context of 
heavy ion collision this implies that MCE effects should be stronger for positively charged particles 
than for negatively charged particles, due to the fact that the created system carries positive net- 
charge. 

Previous work suggests that the asymptotic values for the scaled variance are indeed reached 



rather quickly [15] and above results are certainly applicable to large systems expected to be created 
in relativistic heavy ion collisions. 



VI. REMARKS AND CONCLUSION 

Some concluding remarks seems to be in order. Although it might seem inappropriate to use 
the MCE formulation of a hadron resonance gas model for calculation of multiplicity fluctuations 
in heavy ion collisions, as energy and volume cannot be assumed to be the same in all events, it 
should be stressed that GCE and CE still imply a very particular type of heat (and momentum) 
bath, namely an infinite (and ideal) one. This assumption seems to us even less appropriate. Also 
the MCE is often understood as the ensemble with energy (and charge), however not momentum 
conservation. It is usually assumed that taking momentum conservation into account will not affect 
fluctuations in the large volume limit. We have shown [l^] in a recent paper that this is indeed 
the case, when one assumes information about all produced particles. However for calculations of 
multiplicity fluctuations in arbitrary finite subsystems in momentum space all kinematic conservation 
laws need to be taken into account. 

In a realistic heavy ion experiment it seems impossible to measure the entire final state of each 
collision. The observed subsystem could therefore be seen as effectively embedded into a (possibly 
much larger) heat, charge, and momentum bath. Sometimes it is therefore argued that, when in- 
vestigating only a small part of a statistical system (canonical or micro-canonical), one can ignore 
correlations of the subsystem under investigation with the remaining system. This argument is often 
applied when considering yields and/or fiuctuations in a limited segment of momentum space. More 
precisely, usually the GCE is thought to be the appropriate ensemble to model fiuctuations of particle 
multiplicity or particle ratios found in some mid-rapidity interval In this work we have argued 
that this assumption should be checked carefully. The GCE is only the correct ensemble to choose, 
if heat and charge bath are assumed to be infinite, while the observed subsystem remains finite. 

Based on our previous line of arguments, one would also expect that strong collective longitudinal 
and transverse fiow would lead to a strong correlation of macroscopic subsystems. Longitudinal 
momentum conservation implies that when 'observing' in an event a final state with a certain small 
(large) number of produced particles at very forward rapidity, a similarly small (large) number of 
particles should exist at backward rapidities. Particles in these bins carry substantial longitudinal 
momenta, and hence energy. Modest fiuctuations in their numbers should therefore induce stronger 
fiuctuations in the central rapidity region. The same line of arguments is applicable to the transverse 
momentum dependence. One would therefore expect a similar momentum space dependence of 
experimentally measured charged particle multiplicity fiuctuations as shown in Figs. [2l 



This argument is additionally supported by UrQMD simulations 



291]. In transport calculations 



the produced systems stay far away from global or local equilibrium [30|] and other (dynamical) 



mechanisms might lead to similar effects. On the other hand could one also infer from [29| that 
even in non-equilibrium systems correlations due to exactly enforced conservation laws determine 
the general trend, although transport simulations show, for instance, a very different dependence of 



32| than statistical equilibrium models. This should be 



multiplicity fluctuations on beam energy 31 
subject of further investigation. 

Finally, and most importantly, we want to stress that recently presented preliminary NA49 anal- 
ysis of multiplicity fluctuations in certain rapidity and transverse momentum windows [33] shows 
qualitatively the very same trends as they are suggested by the MCE formulation of the statistical 
model. Data, UrQMD simulations, and the statistical model exhibit suppressed multiplicity fluc- 
tuations when bins with high transverse momentum (or high values of rapidity) are compared to 
bins of same mean multiplicity at lower transverse momentum (or lower values of rapidity). We are 
certainly tempted to interpret this rather unexpected common behavior as a manifestation of energy 
and momentum conservation effects. 



VII. SUMMARY 



We have discussed the effect of momentum space cuts on multiplicity fluctuations in the framework 
of an ideal classical pion gas in the three standard ensembles, GCE, CE, and MCE. Only in the MCE 
we expect a momentum space dependence of multiplicity fluctuations, when comparing intervals of 
same average multiplicity. We have shown that even in the thermodynamic limit energy-momentum 
conservation can leave a sizable effect in the fluctuation pattern. 

In a previous publication we have argued that despite the fact one may expect event-by-event 
fluctuations of the thermal energy, i.e. the part of the total energy which goes into thermal particle 
production rather than collective expansion, these event-by-event fluctuations remain small compared 
to energy fluctuations one would expect from grand canonical and canonical ensembles. In this work 
we have shown that energy and momentum conservation lead to a non-trivial momentum space 
dependence of the fluctuation pattern. This argument seems to be strongly supported by data. 

Above results become all the more interesting when compared to models which seek to describe 
effects beyond our considerations. In fact our calculations suggest a similar strength of respective 
suppression or enhancement as they were predicted as signals for the critical point of strongly in- 
teracting matter, the onset of deconflnement, or generally a possible phase transition. One might 
also be tempted to argue, that enhanced fluctuations around mid-rapidity, when compared to a more 
forward rapidity slice, should be interpreted as a signal of a phase transition from a quark gluon 
plasma to a hadron gas phase, expected to be flrst realized in the presumably hotter and denser 
central rapidity region. However in this case there should be a non-monotonic variation as center 



of mass energy of colliding nuclei is changed. This seems not to be supported by preliminary NA49 
data. 

In summary, above results should be treated as a prediction for general trends of multiplicity 
fluctuations in limited segments of momentum space. The existence of this general behavior should 
be further tested by current experiments. Observation of effects similar to those of Figs. [2] in 
experimental data would, in our opinion, strongly speak in favor of our hypothesis that fluctuations 
of extensive observables are indeed dominated by material and motional conservation laws. 
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APPENDIX A: GLOBALLY CONSERVED QUANTITIES 



Turning now to calculations of cumulants, Eq.( |T8|) . we employ always coordinates most suitable 
to our problem. The invariant phase space element is given by: 



dN 



dN 
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_ _ - ^ exp ( - ] , (Al) 

d^p rriT drriT dy dcj) dpx dy d(f) (27r)'^ \ T J ' 

where the single particle energy e = ttlt coshy, its longitudinal momentum pz = mx sinhy, transverse 

mass mip = Px + transverse momentum pip = pi +Py, and rapidity y = tanh {pz/e). Additionally 

we employ spherical coordinates: 



dN _ .^^ 2 

d^p ^ d(j) d9 dp 



(A2) 



For clarity we consider explicitely a few terms, not given in 17|, here. The total energy density is 
given by the sum over individual contributions of all particle species k: 
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The diagonal energy element ^2 ' is given by: 
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Additionally we define the diagonal momentum correlation terms, with pz = pcos6: 
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Due to spherical symmetry in momentum space we find k^""^"^ ~ '^^^'^^ 



Pz,Pz 



(A5) 



Correlation terms 



of odd order in one of the momenta are identical to zero. As an example we find for correlations 
between energy and longitudinal momentum: 
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since the integral over the polar angle JJ^sin6'cos0 = 0. Similarly we find k 
Additionally kF^ = 0, etc., since for a static source (P) = 0. 
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APPENDIX B: TRANSVERSE MOMENTUM SEGMENT 



The average particle number density of tt in a segment of transverse momentum ApT is given 
by Eq. (fT8!) . i.e. the first derivative of the CGF with respect to (j)NQ = <pNp^ at the origin: 
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Please note that k^*"^ = f^p^dpx dN/dpx = {Np^). Correlations of tx~ in a segment Apj- with 
globally conserved energy are given by double differentiation of \E' with respect to ipNp^ and 0^;, 
thus: 
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Correlations between conserved momenta and particles in Apx, given by the elements k. 
and 1^2'"^ are identical to zero, due to symmetry in azimuthal angle (p for the first two, and due 



to an antisymmetric rapidity integral for the last. Therefore, all elements involving an odd order in 
one of the momenta in Eq. flMl) are equal to zero. The determinant of Eq. flMl) thus factorizes, similar 
to Eq. (l26l) . into a product of {k,2^'^'"Y the determinant of a 3 x 3 sub-matrix involving only 
terms containing Np^,E,Q. Hence momentum conservation drops out when calculating Eq.f l22l) . 
However the strength of correlations between particle number Np^ and globally conserved energy E 
will depend on the position of the segment Apx- Thus using Eqs. (l22i) and ( l24l) . one can express the 
width of the MCE multiplicity distribution fl2Tl) by: 
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In Boltzmann approximation, we find from Eq. ffTSj) . Kg 



K 



Pt 



g/tf^"", where we 



have defined the acceptance q 



N, 



Pt 



/ K^*'' . However, when observing a fraction q of the particle 



density, one does not necessarily observe the same fraction q of the energy density (E^) carried by 
7i~, and thus k^^'^'^ ^ q{EJ). Therefore depending on the location of Ap^, our detector sees a larger 
(smaller) fraction of the total energy, which leads to smaller (larger) particle number fluctuations, see 
Fig. [21, left panel. One can easily verify that setting Kg*"^' = q{EJ) in Eq. flB3p . leads to acceptance 
scaling, Eq.dSD, u^;'^ = l + q {ujZ' - 1). 



APPENDIX C: RAPIDITY SEGMENT 



The average particle number density of vr in a rapidity interval Ay is given by: 



K 



Ny 
1 



d 



dmi 



dy 



dN 



drriT dy 



9 e ^ 
i2ny 



Ay 
2 



T dy exp{^-— cosh {y)j {^—j + 2— cosh ?/ + 2cosh y 



Ay 



(CI) 



Please note that = J^^ dy dN/ dy = (Ny) . Correlations of particles in Ay with globally conserved 
energy are given by: 
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The correlation term of particles in Ay with globally conserved longitudinal momentum Pz reads: 
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Thus the element k^^''^'' in the matrix (12^ is no n- vanishing, and longitudinal momentum (P^) conser- 
vation seems to affects correlations between particles in a segment Ay and the remaining system. In 



contrast to that further elements are equal to zero, k 



K 



Ny,Py 



0, and Px and Py conservation 



have no additional effect. When momentum conservation is taken into account the scaled variance 
can be calculated from Eq. (125]l : 
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Similarly to the previous section, we find a large (small) k^^'^^ leads to small (large) fluctuations, 
see Fig. [21 right panel. When intervals symmetric in rapidity are assumed, e.g. Ay = [—yi,yi\, or 
Ay = [—y2, — ?/i]U[?/i, 7/2], correlations between particle number and momentum disappear, n 



Ny,Pz 



0, 



and Eq. flC4l) reduces to Eq. flBSp . and momentum conservation does not play a role. Equally when 
disregarding longitudinal momentum conservation the same arguments as those of Appendix [B] apply 
and Eq.( ]B3[ ) holds, however the effect is much weaker, see Fig. |3l 



APPENDIX D: AZIMUTHAL ANGLE SEGMENT 



Th average particle number in Ac/), while integrating over all pt and y is simply a fraction q 



A0/27r of the total yield (N^n)- Therefore k^''' = qn^^'' . Equally, the energy carried by vr in this 



E,N^ 



interval is k^'"'* = q{E^). Due to symmetry around y = 0, we find additionally 1^2'*"^'' = 0- However 
for the transverse momenta Px = Pr cos 0, and Py = pt sin the correlation with A'^^ is generally 
no n- zero. 



Nj„Px 



d 



)COS( 



d 



-'Px 



9 



<t>,=o 



(27r 



2g e 



A(j> 



(2vr) 



rn^ T mT Kr^/o ( ^ 



00 00 

dpT 
-00 



dy P:c 



dN 



5/2 



VT 



-1 



sm I 



dpT dy 
(Pt) ■ 



A<j> 



(Dl) 



Similarly we find k^*^'^" = — (27r) ""^[cos^]^^ (pr)- Unlike in the previous sections there is no 
particular dependence of the position of the interval A0. However in general there is a dependence. 
When momentum conservation is taken into account Eq. fl22l) can be calculated from Eg. (1251) : 
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where we have used k^"'^"" 



^Pv^Py 



. As mentioned before there is no particular dependence of Kg 
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and K^'^'^ on the position of A0. However we have a term ^/t^'*'^'' j + 
a continuous interval = 02] this terms reads: 
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This term is evidently positive, hence fluctuations are suppressed. One can easily verify that when 
one takes A^^ = ^2] U [cpi + vr, 02 + vr], i.e. two opposite slices in azimuthal angle, the correlation 



disappears, k 



N4„Px 



N^,Py 



0, and one returns to acceptance scaling, Eq.( 



[1] J. Cleymans, H. Oeschler, K. Redlich and S. Wheaton, Phys. Rev. C 73 (2006) 034905. 
[2] F. Becattini, J. Manninen and M. Gazdzicki, Phys. Rev. C 73 (2006) 044905. 
[3] A. Andronic, P. Braun-Munzinger and J. Stachel, Nucl. Phys. A 772 (2006) 167. 
[4] J. Letessier and J. Rafelski, |arXiv:nucl-th /0504028 1 

[5] B. Lungwitz et al. [NA49 Collaboration], PoS C FRNC2006 (2006) 024. 

[6] V. V. Begun, M. Gazdzicki, M. I. Gorenstein, M. Hauer, V. P. Konchakovski and B. Lungwitz, Phys. 

Rev. C 76 (2007) 024902. 
[7] V. V. Begun, M. Gazdzicki, M. I. Gorenstein and O. S. Zozulya, Phys. Rev. C 70, 034901 (2004). 
[8] H. Heiselberg, Phys. Rept. 351 (2001) 161; S. Jeon and V. Koch, in Quark-Gluon Plasma 3, edited by 

R.C. Hwa and X.-N. Wang (World Scientific, Singapore, 2004), p.430. 
[9] M. Gazdzicki, M. I. Gorenstein and S. Mrowczynski, Phys. Lett. B 585 (2004) 115; M. I. Gorenstein, 

M. Gazdzicki and O. S. Zozulya, Phys. Lett. B 585 (2004) 237. 
[10] I. N. Mishustin, Phys. Rev. Lett. 82 (1999) 4779; H. Heiselberg and A. D. Jackson, Phys. Rev. C 63 

(2001) 064904. 

[11] M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81 (1998) 4816; Phys. Rev. D 60 
(1999) 114028; M. Stephanov, Acta Phys. Polon. B 35 (2004) 2939. 



[12] S. Haussler, S. Scherer and M. Bleicher, [arXiv:hep-ph/0702188[ 



[13] T. K. Navak. larXiv:0706:2708^ [nucl-ex], 

[14] V. V. Begun, M. I. Gorenstein, M. Hauer, V. P. Konchakovski and O. S. Zozulya, Phys. Rev. C 74 
(2006) 044903. 

[15] F. Becattini, A. Keranen, L. Ferroni and T. Gabbriellini, Phys. Rev. C 72 (2005) 064904. 
[16] M. I. Gorenstein, M. Hauer and D. O. Nikolajenko, Phys. Rev. C 76, 024901 (2007). 



[17] M. Hauer, V. V. Begun and M. I. Gorenstein, arXiv:0706.3290 [nucl-th]. 

[18] V. V. Begun and M. I. Gorenstein, Phys. Rev. C 73 (2006) 054904; V. V. Begun, M. I. Gorenstein, 

A. P. Kostyuk and O. S. Zozulya, J. Phys. G 32 (2006) 935. 
[19] K. Redhch and L. Turko, Z. Phys. C 5 (1980) 201; L. Turko, Phys. Lett. B 104 (1981) 153; R. Hagedorn 

and K. Redhch, Z. Phys. C 27 (1985) 541. 
[20] J. Cleymans, K. Redhch and E. Suhonen, Z. Phys. C 51 (1991) 137; A. Keranen and F. Becattini, 

Phys. Rev. C 65 (2002) 044901; F. Becattini and U. W. Heinz, Z. Phys. C 76, 269 (1997). 
[21] O. Fochler, S. Vogel, M. Bleicher, C. Greiner, P. Koch-Steinheimer and Z. Xu, Phys. Rev. C 74, 034902 

(2006). 

[22] C. M. Ko, V. Koch, Z. w. Lin, K. Redhch, M. A. Stephanov and X. N. Wang, Phys. Rev. Lett. 86, 
5438 (2001). 

[23] S. Jeon, V. Koch, K. Redhch and X. N. Wang, Nucl. Phys. A 697, 546 (2002). 

[24] F. Becattini and L. Ferroni, Eur. Phys. J. C 35 (2004) 243; Eur. Phys. J. C 38 (2004) 225. 

[25] J. Letessier, A. Tounsi and J. Rafelski, Phys. Lett. B 475 (2000) 213; J. Rafelski and J. Letessier, Phys. 

Rev. Lett. 85 (2000) 4695. 
[26] P. Koch, B. Muher and J. Rafelski, Phys. Rept. 142 (1986) 167. 
[27] V. V. Begun, M. I. Gorenstein and O. S. Zozulya, Phys. Rev. C 72 (2005) 014902. 
[28] S. Jeon and V. Koch, Phys. Rev. Lett. 83, 5435 (1999); G. Torrieri, S. Jeon and J. Rafelski, Phys. Rev. 

C 74 (2006) 024901; G. Torrieri, [a?Xiv:nucl-th/0702062| 



[29] B. Lungwitz and M. Bleicher, arXiv:0707.r788] [nucl-th]. 

[30] E. L. Bratkovskaya, W. Gassing, C. Greiner, M. Effenberger, U. Mosel and A. Sibirtsev, Nucl. Phys. A 

681 (2001) 84; L. V. Bravina et al, Phys. Rev. C 60, 024904 (1999). 
[31] V. P. Konchakovski, M. 1. Gorenstein and E. L. Bratkovskaya, Phys. Lett. B 651 (2007) 114. 
[32] V. P. Konchakovski, M. I. Gorenstein and E. L. Bratkovskava. ll arXiv:0704.1831l [nucl-th]. 
[33] B. Lungwitz et al. [NA49 Collaboration], |ax2£iYiQ709J646J [nucl-ex]. 



